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C^l ABSTRACT 

Q ■ We combine the 2MASS extended source catalogue and the 2dF galaxy redshift sur- 

?-) | vey to produce an infrared-selected galaxy catalogue with 17,173 measured redshifts. 

c/3 , We use this extensive dataset to estimate the galaxy luminosity functions in the J- 

and Kg-bands. The luminosity functions are fairly well fit by Schechter functions with 
parameters M* - 5 log h = -22.36 ± 0.02, aj = -0.93 ± 0.04, $* = 0.0104 ± 0.0016ft 3 
Mpc" 3 in the J-band and Af* g - 5 log ft = -23.44 ± 0.03, a Ks = -0.96 ± 0.05, 
$k g = 0.0108 ± 0.0016ft 3 Mpc~ 3 in the K s -band (2MASS Kron magnitudes). These 

^ . parameters are derived assuming a cosmological model with ilg = 0.3 and Ao = 0.7. 

With datasets of this size, systematic rather than random errors are the dominant 
source of uncertainty in the determination of the luminosity function. We carry out a 
careful investigation of possible systematic effects in our data. The surface brightness 
distribution of the sample shows no evidence that significant numbers of low surface 
brightness or compact galaxies are missed by the survey. We estimate the present-day 
distributions of bj— Ks and J— Kg colours as a function of absolute magnitude and 
use models of the galaxy stellar populations, constrained by the observed optical and 
infrared colours, to infer the galaxy stellar mass function. Integrated over all galaxy 
masses, this yields a total mass fraction in stars (in units of the critical mass density) 
of fWsft = (1.6 ±0.24) x 10" 3 for a Kennicutt IMF and O stars /i = (2.9 ±0.43) x 10~ 3 
for a Salpeter IMF. These values are consistent with those inferred from observa- 
tional estimates of the total star formation history of the universe provided that dust 
extinction corrections are modest. 

Key words: galaxies: luminosity function estimators 



* This publication makes use of data products from the Two Mi- Center / California Institute of Technology, funded by the National 

cron All Sky Survey (2MASS), which is a joint project of the Uni- Aeronautics and Space Administration and the National Science 

versity of Massachusetts and the Infrared Processing and Analysis Foundation. 
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1 INTRODUCTION 

The near-infrared galaxy luminosity function is an impor- 
tant characteristic of the local galaxy population. It is a 
much better tracer of evolved stars, and hence of the total 
stellar content of galaxies, than optical luminosity functions 
which can be dominated by young stellar populations and 
are also strongly affected by dust extinction. Hence, infrared 
luminosities can be much more directly related to the un- 
derlying stellar mass of galaxies and so knowledge of the 
present form and evolution of the infrared galaxy luminos- 
ity function places strong constraints on the history of star 
formation in the universe and on galaxy formation models 
(e.g. Cole et al. 2000 and references therein). 

The local K-band luminosity function has been es- 
timated from optically selected samples by Mobasher, 
Sharpies and Ellis (1993), Szokoly et al. (1998) and Loveday 
(2000) and from K-band surveys by Glazebrook et al. (1995), 
and Gardner et al. (1997). The existing K-band surveys are 
small. The largest, by Gardner et al., covers only 4 deg 2 and 
contains only 510 galaxies. The recent survey of Loveday 
covers a much larger solid angle. In this survey the redshifts 
were known in advance of measuring the K-band magnitudes 
and this was exploited by targetting bright and faint galax- 
ies resulting in an effective sample size much larger than the 
345 galaxies actually measured. However, like all optically 
selected samples, it suffers from the potential problem that 
galaxies with extremely red infrared to optical colours could 
be missed. In this paper we combine the 2-Micron All Sky 
Survey (2MASS) with the 2dF galaxy redshift survey (2dF- 
GRS) to create an infrared selected redshift survey subtend- 
ing 2151.6 deg 2 . Currently the sky coverage of both surveys 
is incomplete, but already the overlap has an effective area 
of 619 deg 2 . Within this area the redshift survey is com- 
plete to the magnitude limit of the 2MASS catalogue and so 
constitutes a complete survey which is 50 times larger than 
the previous largest published infrared selected redshift sur- 
vey. A new catalogue of a similarly large area, also based 
on 2MASS, has very recently been analysed by Kochanek et 
al. (2001). They adopt isophotal rather than total magni- 
tudes and concentrate on the dependence of the luminosity 
function on galaxy morphology. 

In Section 2.1 we briefly describe the relevant proper- 
ties of the 2dFGRS and 2MASS catalogues. Section 2.2 is 
a detailed examination of the degree to which the matched 
2MASS-2dFGRS galaxies are a complete and representative 
subset of the 2MASS catalogue. Section 2.3 examines the 
calibration of the 2MASS total magnitudes and Section 2.4 
demonstrates that the 2MASS catalogue and the inferred 
luminosity functions are not affected by surface brightness 
selection effects. In Section 3 we present the method by 
which we compute k-corrections and evolutionary correc- 
tions and relate the observed luminosities to the underlying 
stellar mass. The estimation methods and normalization of 
the luminosity functions are described briefly in Section 4. 
Our main results are presented and discussed in Section 5. 
These include estimates of the J and Ks (K-short) lumi- 
nosity functions, the bj— Ks and J— Ks colour distributions 
as a function of absolute magnitude and the distribution 
of spectral type. We also estimate the stellar mass function 
of galaxies, which can be integrated to infer the fraction of 



baryons in the universe which are in the form of stars. We 
conclude in Section 6. 



2 THE DATASET 

The data that we analyze are the extended source cata- 
logue from the second incremental release of the 2-Micron 
All Sky Survey (2MASS http://pegasus.phast.umass.edu) 
and the galaxy catalogue of the 2dF galaxy redshift survey 
(2dFGRS http://www.mso.anu.edu.au/2dFGRS). Here, we 
present the relevant properties of these two catalogues and 
investigate their selection characteristics and level of com- 
pleteness. 

2.1 Selection Criteria 

The 2MASS is a ground-based, all-sky imaging survey in the 
J, H and Ks bands. Details of how extended sources are iden- 
tified and their photometric properties measured are given 
by Jarrett et al. (2000). The detection sensitivity (10cr) for 
extended sources is quoted as 14.7, 13.9 and 13.1 magnitudes 
in J, H and Ks respectively. The complete survey is expected 
to contain 1 million galaxies of which approximately 580,000 
are contained in the second incremental data release made 
public in March 2000. 

The 2dFGRS is selected in the photographic bj 
band from the APM galaxy survey (Maddox et al. 
1990a, 1990b, 1996) and subsequent extensions to it, that 
include a region in the northern galactic cap (Maddox 
et al. in preparation). The survey covers approximately 
2151.6 deg 2 consisting of two broad declination strips. 
The larger is centred on the SGP and approximately cov- 
ers -22°.5><5>-37°5, 21 h 40 m <a<3 h 30 m ; the smaller strip 
is in the northern galactic cap and covers 2°.5><5> — 7°.5, 
9 h 50 m <a<14 h 50 m . In addition, there are a number of ran- 
domly located circular 2-degree fields scattered across the 
full extent of the low extinction regions of the southern APM 
galaxy survey. There are some gaps in the 2dFGRS sky cov- 
erage within these boundaries due to small regions that have 
been excluded around bright stars and satellite trails. The 
2dFGRS aims to measure the redshifts of all the galaxies 
within these boundaries with extinction-corrected bj mag- 
nitudes brighter than 19.45. When complete, at the end of 
2001, 250,000 galaxy redshifts will have been measured. In 
this paper we use the 140,000 redshifts obtained prior to 
September 2000. 

The overlap of the two surveys is very good. There are 
some gaps in the sky coverage due to strips of the sky that 
were not included in the 2MASS second incremental release, 
but overall a substantial fraction of the 2151.6 deg 2 of the 
2dFGRS is covered by 2MASS. The homogeneity and exten- 
sive sky coverage of the combined dataset make it ideal for 
studies of the statistical properties of the galaxy population. 

2.2 The Completeness of the Matched 
2MASS-2dFGRS Catalogue 

Here we consider whether all the 2MASS galaxies within 
the 2dFGRS survey region have 2dFGRS counterparts and 
assess the extent to which the fraction of galaxies with mea- 
sured redshifts represents an unbiased sub-sample. 
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Figure 1. The distribution of angular separation, 0, for matched 
2MASS-2dFGRS galaxies. The solid histogram is the distribution 
for the whole catalogue and the dotted histogram for the subset 
of 2MASS galaxies with semi-major axes larger than 12 arcsec. 



The astrometry in both 2MASS and 2dFGRS is, in gen- 
eral, very good and it is an easy matter to match objects in 
the two catalogues. We choose to find the closest pairs within 
a search radius equal to three quarters of the semi-major axis 
of the J-band image (denoted j_r_e in the 2MASS database). 
Scaling the search radius in this way helps with the matching 
of large extended objects. This procedure results in the iden- 
tification of 2dFGRS counterparts for 40,121 of the 2MASS 
objects, when at random one would only expect to find a 
handful of such close pairs. Moreover, the distribution of 
separations shown in Fig. 1 peaks at 0.5 arcsec, with only 
3% having separations greater than 3 arcsec. A significant 
part of this tail comes from the most extended objects as 
is evident from the dotted histogram in Fig. 1 which shows 
objects with semi-major axes larger than 12 arcsec. Thus, 
we can be very confident in these identifications. 

The 40,121 2MASS objects for which we have found se- 
cure 2dFGRS counterparts amount to 88.6% of the 2MASS 
extended sources that fall within the boundary of the 2dF- 
GRS. As discussed below, a more restrictive criterion that 
includes only sources fainter than J=12 that are confidently 
classified as galaxies by 2MASS, increases the fraction with 
2dFGRS matches to 90.7%. The remaining 9.3% are missed 
for well understood reasons (star-galaxy classification: 4.6%; 
merged or close images: 4.4%; miscellaneous: 0.27%), none 
of which ought to introduce a bias. This is confirmed explic- 
itly, in the middle row of Fig. 2, by the close correspondence 
between the photometric properties of the missed 9.3% and 
those of the larger matched sample. Hence, in estimating 
luminosity functions no significant bias will be introduced 
by assuming the matched sample to be representative of the 
full population. Furthermore, the distribution shown in the 
bottom row of Fig. 2 shows that the subset of 17,173 galaxies 
for which we have measured redshifts is a random sample of 
the full matched 2MASS-2dFGRS catalogue. This summary 
is the result of a thorough investigation, which we describe 
in the remainder of this section, into the reasons why 11.4% 



of the 2MASS sources are missed and whether their omission 
introduces a bias in the properties of the matched sample. 

We first consider objects in the 2MASS catalogue which 
based on their images and colours are not confidently clas- 
sified as galaxies. In the 2MASS database a high e_score or 
g_score indicates a high probability that the object is either 
not an extended source or not a galaxy. A ccJlag/O indicates 
an artifact or contaminated and/or confused source. For de- 
tailed definitions of these parameters we refer the reader to 
Jarrett et al. (2000). Rejecting all objects which have either 
e_score>1.4, g_score>1.4 or cc_flag/0 removes just 6.7% of 
the total. However, removing these reduces significantly the 
fraction of the 2MASS sample that does not match with 
the 2dFGRS catalogue, from 11.4% to 9.6%. Thus, it is 
likely that about 30% of the 2MASS objects which have 
e_score>1.4, g_score>1.4 or ccJlag/O are not galaxies. 

The 2MASS may contain a tail of very red objects 
that are too faint in the bj-band to be included in the 
b.j< 19.45 2dFGRS sample. Fig. 3 shows the distribution of 
bj-J colours for the matched objects with J<14.7. (Here, 
the J-band magnitude we are using is the default magni- 
tude denoted j_m in the 2MASS database. In Section 2.3 
we will consider the issue of what magnitude definition is 
most appropriate for estimating the luminosity function.) 
The vertical dashed line indicates the colour at which this 
sample starts to become incomplete due to the bj< 19.45 
magnitude limit of 2dFGRS. The colour distribution cuts 
off sharply well before this limit, suggesting that any tail of 
missed very red objects is extremely small. In other words 
the 2dFGRS is sufficiently deep that even the reddest ob- 
jects detected at the faintest limits of 2MASS ought to be 
detected in 2dFGRS. 

In the top row of Fig. 2 we compare the distributions 
of magnitude, colour and surface brightness for the matched 
and missed 2MASS objects. In general, the properties of the 
missed subset overlap well with those of the much larger 
matched subset. However, we do see that the distributions 
for missed objects contain tails of bright and blue objects. It 
is quite likely that this is due to the 2MASS extended source 
catalogue being contaminated by a small population of sat- 
urated or multiple stars. The dotted histograms in the top 
row of Fig. 2 show the distributions of magnitude, colour 
and surface brightness for the bright subset of the missed 
objects with J<12. Here we clearly see bimodal colour and 
surface brightness distributions. The blue peak of the colour 
distribution is consistent with that expected for stars (see 
Jarrett et al. 2000). Excluding these bright, J<12, objects 
which are clearly contaminated by stars reduces the fraction 
of missed 2MASS objects from 9.6% to 9.3%. The magni- 
tude, colour and surface brightness distributions for this re- 
maining 9.3% are shown in the middle panel of Fig. 2. We 
see that the missed objects are slightly under-represented at 
the faintest magnitudes and also slightly bluer on average 
than the matched sample, while the distribution of surface 
brightness is almost indistinguishable for the two sets of ob- 
jects. These differences are small and so will introduce no 
significant bias in our luminosity function estimates. 

To elucidate the reasons for the remaining missed 9.3% 
of 2MASS objects we downloaded 100 lxl arcmin images 
from the STScI Digitized Sky Survey (DSS) centred on the 
positions of a random sample of the missed 2MASS objects. 
In each image we plotted a symbol to indicate the position of 
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Figure 2. The distribution of J-band apparent magnitude, J— Kg colour and J-band surface brightness, juj, for various sub-samples of 
the 2MASS catalogue. Here, the measure of surface brightness used is simply /Ltj = J — 51og 10 r, where J is the Kron magnitude and r 
the Kron semi-major axis in arc seconds (j_m_e and j_r_e in the 2MASS database). In all three rows, the thick solid histograms are the 
distributions for 2MASS objects that are matched with 2dFGRS galaxies. The light solid histograms in the top row are the 11.4% of 
2MASS galaxies that are not matched with 2dFGRS galaxies. Poisson errorbars are shown on these histograms. The dashed histograms 
are for the bright sub-sample with J<12. In the middle row, the light histograms show the distributions for the 9.3% of 2MASS galaxies 
fainter than J=12 and satisfying the additional image classification constraints discussed in the text that are not matched with 2dFGRS 
galaxies. In the bottom row, the light histograms show the distributions of the 42.8% of the matched 2MASS-2dFGRS galaxies for which 
redshifts have been measured. The values in each histogram are the fraction of the corresponding sample that falls in each bin. 
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Figure 3. The solid histogram shows the distribution of bj-J 
colours for 2MASS galaxies selected to have J<14.7. (Here, we 
use the 2MASS default magnitude, denoted j_m in the 2MASS 
database.) The vertical dashed line indicates the colour at which 
this sample starts to become incomplete due to the bj< 19.45 
magnitude limit of 2dFGRS. 



any 2dFGRS galaxies within the lxl arcmin field. We also 
plotted symbols to indicate the positions and classifications 
of all images identified in the APM scans from which the 
2dFGRS catalogue was drawn, down to a magnitude limit 
of bj~ 20.5. These images are classified as galaxies, stars, 
merged images (galaxy+galaxy, galaxy+star or star+star) 
or noise. This set of plots allows us to perform a census of 
the reasons why some 2MASS objects are not present in the 
2dFGRS survey. 

The main cause for the absence of 2MASS objects in 
the 2dFGRS is that the APM has classified these objects as 
stars. These amount to 49.5% of the missed sample (4.6% 
of the full 2MASS sample). In some cases, the DSS image 
shows clearly that these are stars and in others that they 
are galaxies. However, the majority of these objects cannot 
easily be classified from the DSS images. Thus, they could 
be galaxies that the APM has falsely classified as stars or 
stars that 2MASS has falsely classified as galaxies. The first 
possibility is not unexpected since the parameters used in 
the APM star-galaxy separation algorithm were chosen as 
a compromise between high completeness and low contami- 
nation such that the expected completeness is around 95% 
with 5% stellar contamination (Maddox et al. 1990a). It is 
hard to rule out the possibility that this class of object does 
not include a substantial fraction of stars, but if so, their 
presence appears not to distort the distribution of colours 
shown in Fig 3. Another 47.6% of the random sample (4.4% 
of the full 2MASS sample) are classified by the APM as 
mergers or else consist of two close images in the DSS but 
are classified by the APM as a single galaxy offset from the 
2MASS position. The remaining 2.9% of the random sample 
(0.27% of the full 2MASS sample) are missed for a variety of 
reasons including proximity to the diffraction spikes of very 
bright stars and poor astrometry caused by the presence of 
a neighbouring unclassified image. 
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Figure 4. A comparison of the 2MASS default, Kron and ex- 
trapolated magnitudes in the J and Kg bands. The dots are the 
measured values for each of the galaxies in the matched 2MASS- 
2dFGRS catalogue. The solid and dotted lines indicate the me- 
dian, 10 and 90 percentiles of the distribution. 



2.3 2MASS Magnitude Definitions and 
Calibration 

The 2MASS extended source database provides a large selec- 
tion of different magnitude measurements. In the previous 
section we used the default magnitudes (denoted j_m and 
k_m in the 2MASS database). These are magnitudes de- 
fined within the same circular aperture in each waveband. 
For galaxies brighter than Ks=14, the aperture is the cir- 
cular Ks-band isophote of 20 mag arcsec~ 2 and for galax- 
ies fainter than Ks=14 it is the circular J-band isophote of 
21 mag arcsec -2 . These are not the most useful definitions 
of magnitude for determining the galaxy luminosity func- 
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Figure 5. Comparison of 2MASS Kron and extrapolated magnitudes with the independent measurements of Loveday (2000). The left 
hand panels are for the Kg-band Kron and extrapolated magnitudes (k_m_e and k_m_ext in the 2MASS database). The right hand panels 
show Kron and extrapolated magnitudes inferred from the 2MASS J-band Kron and extrapolated magnitudes and the measured default 
aperture J— Kg colours (j_m_e— j_m+k_m and j_m_ext— j_m+k_m in the 2MASS database variables). The horizontal errorbars show the 
measurement errors quoted by Loveday (2000). The solid lines show simple least squares fits. The slopes and zero-point offsets of these fits 
and the rms residuals about the fits are indicated on each panel. The inset plots show the distribution of residual magnitude differences. 



tion. Since we are interested in measuring the total lumi- 
nosity and ultimately the total stellar mass of each galaxy, 
we require a magnitude definition that better represents the 
total flux emitted by each galaxy. We consider Kron magni- 
tudes (Kron 1980) and extrapolated magnitudes. Kron mag- 
nitudes (denoted j_m_e and k_m_e in the 2MASS database) 
are measured within an aperture, the Kron radius, defined 
as 2.5 times the intensity-weighted radius of the image. The 
extrapolated magnitudes (denoted j_m_ext and k_m_ext in 
the 2M ASS database) are defined by first fitting a modified 
exponential profile, f(r) = fo exp[— (ar) 1 ^], to the image 
from 10 arcsec to the 20 mag/arcsec 2 isophotal radius, and 



extrapolating this from the Kron radius to 4 times this ra- 
dius or 80 arcsec if this is smaller (Jarrett private commu- 
nication). Note that improvements are being made to the 
extended source photometry algorithms developed and em- 
ployed 2MASS team and so in the final 2MASS data release 
the definitions of the Kron and extrapolated magnitudes 
may be slightly modified (Jarrett private communication). 

Fig. 4 compares the default, Kron and extrapolated 
magnitudes in the J and Ks bands for the matched 2MASS- 
2dFGRS catalogue. The upper panel shows that while the 
median offset between the J-band isophotal default magni- 
tudes and the pseudo-total Kron magnitudes is small there is 
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a large spread with some galaxies having default magnitudes 
more than 0.5 magnitudes fainter than the Kron magnitude. 
The Kron magnitudes are systematically fainter than the 
extrapolated magnitudes by between approximately 0.1 and 
0.3 magnitudes. This offset is rather larger than expected: if 
the Kron radius is computed using a faint isophote to define 
the extent of the image from which the intensity weighted 
radius is measured, then the Kron magnitudes should be 
very close to total. For an exponential light profile (f3 = 1), 
the Kron radius should capture 96% of the flux, while for 
an r 1 / 4 law (/3 = 4), 90% of the flux should be enclosed. 
In other words, the Kron magnitude should differ from the 
total magnitude by only 0.044 and 0.11 magnitudes in these 
two cases. However, the choice of isophote is a compro- 
mise between depth and statistical robustness. In the case 
of the 2MASS second incremental release, an isophote of 
21.7(20.0) mag arcsec -2 in J(Ks) was adopted (Jarrett pri- 
vate communication). These relatively bright isophotes, par- 
ticularly the Ks-band isophote, could lead to underestimates 
of the Kron radii and fluxes for lower surface brightness ob- 
jects and plausibly accounts for much of the median offset 
of 0.3 magnitudes seen in Fig. 4 between the Ks-band Kron 
and extrapolated magnitudes. This line of reasoning favours 
adopting the extrapolated magnitudes as the best estimate 
of the total magnitudes, but, on the other hand, the ex- 
trapolated magnitudes are model-dependent and have larger 
measurement errors. 

To understand better the offset and scatter in the 
2MASS magnitudes we have compared a subset of the 
2MASS data with the independent K-band photometry of 
Loveday (2000) . The pointed observations of Loveday have 
better resolution than the 2MASS images and good signal- 
to-noise to a much deeper isophote. This enables accurate, 
unbiased Kron magnitudes to be measured. Note that the 
offset between the 2MASS Ks-band and the standard K- 
band used by Loveday is expected to be almost completely 
negligible (see Carpenter 2001). The left hand panels of 
Fig. 5 compare these measurements with the corresponding 
2MASS Kron and extrapolated magnitudes. The right hand 
panels show Ks-band Kron and extrapolated magnitudes 
computed by taking the 2MASS J-band Kron and extrap- 
olated magnitudes and subtracting the J— Ks colour mea- 
sured within the default aperture. These indirect estimates 
are interesting to consider as they combine the profile infor- 
mation from the deeper J-band image with the J— Ks colour 
measured within the largest aperture in which there is good 
signal-to-noise. The straight lines plotted in Fig. 5 show sim- 
ple least squares fits and the slope and zero-point offset of 
these fits are indicated on each panel along with bootstrap 
error estimates. Also shown in the inset panels is the distri- 
bution of residual magnitude differences about each of the 
fits and a gaussian fit to this distribution. The rms of these 
residuals and a bootstrap error estimate is also given in each 
panel. 

From these comparisons we first see that all the fits 
have slopes entirely consistent with unity, but that their ze- 
ropoints and scatters vary. The zero-point offsets, A^ r ° n , be- 
tween both the 2MASS Kron magnitude measurements and 
those of Loveday confirm that the 2MASS Kron magnitudes 
systematically underestimate the galaxy luminosities. In the 
case of the direct Ks-band 2MASS magnitudes the offset is 
A^ r , on = 0.164 magnitudes. In the case of the Kron magni- 



tudes inferred from the deeper J-band image profiles, the off- 
set is reduced to A^j on = 0.061 magnitudes. Conversely the 
2MASS extrapolated magnitudes are systematically brighter 
than the Loveday Kron magnitudes by — A°*j trap ' = 0.137 
and 0.158 magnitudes, where one would expect an offset of 
only Aktoii = 0.044 to 0.11 due to the difference in definition 
between ideal Kron and true total magnitudes. For both es- 
timates of the extrapolated magnitude and for the directly 
estimated Kron magnitude the scatter about the correlation 
is approximately 0.14 magnitudes and we note a slight ten- 
dency for the scatter to increase at faint magnitudes. The 
magnitude estimate that best correlates with the Loveday 
measurements is the Kron magnitude estimated from the 
2MASS J-band Kron magnitude and the default aperture 
J— Ks colour. Here the distribution of residuals has a much 
reduced scatter of only 0.1 magnitudes and has very few 
outliers. 

Our conclusions from the comparison of Kron magni- 
tudes is that it is preferable to adopt the Ks-band magni- 
tude inferred from the J-band Kron or extrapolated mag- 
nitude by converting to the Ks-band using default aper- 
ture colour, rather than to use the noisier and more bi- 
ased direct Ks-band estimates. With this definition, we find 
that the 2MASS Kron magnitudes slightly underestimate 
the galaxy luminosities while the extrapolated magnitudes 
slightly overestimate the luminosities, particularly at faint 
fluxes. We will present results for both magnitude defini- 
tions, but we note that to convert to total magnitudes we es- 
timate that the 2MASS Kron magnitudes should be bright- 
ened by A^° n + A Kr on = 0.1-0.17 magnitudes and the 
extrapolated magnitudes dimmed by — A™j trap ' — AKron = 
0.05-0.11 magnitudes. 

2.4 Completeness of the 2MASS Catalogue 

Here we define the magnitude limited samples which we 
will analyze in Section 4 and test them for possible incom- 
pleteness in both magnitude and surface brightness. For the 
Kron and extrapolated magnitudes, the 2MASS catalogue 
has high completeness to the nominal limits of J<14.7 and 
Ks<13.9. However, to ensure very high completeness and 
avoid any bias in our luminosity function estimates, we made 
the following more conservative cuts. For the Kron magni- 
tudes, we limited our sample to either J<14.45 or Ks<13.2, 
and for the extrapolated magnitudes to either J<14.15 or 
Ks<12.9. These choices are motivated by plots such as the 
top panel of Fig. 4. Here the isophotal default magnitude 
limit of J<14.7 is responsible for the right hand edge to the 
distribution of data points. One sees that this limit begins 
to remove objects from the distribution of Kron magnitudes 
for JJJ14.5. An indication that the survey is complete to our 
adopted limits is given by the number counts shown in Fig 6, 
which only begin to roll over at fainter magnitudes. 

More rigorously, we have verified that the samples are 
complete to these limits by examining their V(zi)/V(zm&x,i) 
distributions. Here, Zi is the redshift of a galaxy in the sam- 
ple, Zmax.i is the maximum redshift at which this galaxy 
would satisfy the sample selection criteria, and V(z) is the 
survey volume that lies at redshift less than z. If the sam- 
ple is complete and of uniform density, V(zi)/V(2 ma x,») is 
uniformly distributed within the interval to 1. To evaluate 
Zmax we made use of the default k+e corrections described 



© 0000 RAS, MNRAS 000, 000-000 



8 S. Cole et al. 



-0.5 



-1 



-1.5 



'ht - ♦ * ' 




10 



12 



14 



16 



1.5 



-1 



O Gardner et a] (199f 
□ Glazebrook etal (19 




Figure 6. Differential galaxy number counts in the J and Ks 
bands, all with Poisson crrorbars and with a Euclidean slope sub- 
tracted so as to expand the scale of the ordinate. The J and Ks 
counts linked by the solid line are the 2MASS 7 arcsec aperture 
counts of Jarrett et al. (in preparation). The counts linked by 
the dashed and dotted lines are those of the 2MASS-2dFGRS 
redshift catalogue for Kron and extrapolated magnitudes respec- 
tively. The Ks-band magnitudes are those inferred from the J- 
band magnitudes and aperture colours. In the Ks-band these are 
compared with the counts of Gardner et al. (1996) and Glaze- 
brook et al. (1994) as indicated in the figure legend. 



in the following section, but the results are not sensitive to 
reasonable variations in the assumed corrections or in the 
cosmology. The solid histograms in Fig. 7 show these dis- 
tributions for each of our four magnitude limited samples. 
Note that the Ks-band magnitudes are those inferred from 
the J-band magnitudes and aperture colours. The dashed 
histograms in the lower panels show the corresponding dis- 
tributions for the directly measured Ks-band magnitudes. 
In all these cases we have computed V max simply from the 
imposed apparent magnitude limits and have ignored any 
possible dependence of the catalogue completeness on sur- 
face brightness. 

If the samples were incomplete the symptom one would 
expect to see is a deficit in the V/Vmax distributions at 
large V/Vmax and hence a mean (V/Vmax)<0.5. There is 
no evidence for such a deficit in these distributions. In fact 
each has a mean (V/V m ax) slightly greater than 0.5. The 
slight gradient in the V/V ma x distribution is directly re- 
lated to the galaxy number counts shown in Fig 6, which 
are slightly steeper than expected for a homogeneous, non- 
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Figure 7. The distributions of V/V max for our magnitude limited 
samples. The solid histograms in the four panels show the V/ V ma x 
distributions for our J and Ks Kron and extrapolated magnitude 
limited samples. The mean values of (V/V m ax) are indicated on 
each panel. The Ks-band magnitudes are those inferred from the 
J-band magnitudes and aperture colours. The distributions for 
the directly measured Ks-band Kron and extrapolated magni- 
tudes are shown by the dashed histograms in the lower panels. 
The dotted histogram in the top-left panel shows the V/V m ax 
distribution we obtain when attempting to take account of the 
2MASS isophotal diameter and isophotal magnitude limits in es- 
timating the V m ax values. 



evolving galaxy distribution. A similar result has been found 
in the bright bj-band counts (Maddox et al. 1990c). The 
bj-band result has variously been interpreted as evidence 
for rapid evolution, systematic errors in the magnitude cal- 
ibration, or a local hole or underdensity in the galaxy dis- 
tribution (Maddox et al. 1990c; Metcalfe, Fong & Shanks 
1995; Shanks 1990). Here we note that the gradient in the 
^/Knax distributions (and also in the galaxy counts) be- 
comes steeper both as one switches from Kron to the less re- 
liable extrapolated magnitudes and as one switches from the 
J-band data to the lower signal-to-noise Ks-band data. This 
gives strong support to our decision to adopt the Ks-band 
magnitudes derived from the J-band Kron and extrapolated 
magnitudes and aperture J— Ks colours. It also cautions that 
the mean (V/Vmax) >0. 5 cannot necessarily be taken as a 
sign of evolution or a local underdensity, but may instead be 
related to the accuracy of the magnitude measurements. The 
comparison to the observations of Loveday (2000) shows no 
evidence for systematic errors in the magnitudes, but does 
not constrain the possibility that the distribution of magni- 
tude measurement errors may become broader or skewed at 
fainter magnitudes. Such variations would affect the V/Vmax 
distributions and could produce the observed behaviour. We 
conclude by noting that while the shift in the mean ( V/ V max ) 
is statistically significant, it is nevertheless quite small for 
the samples we analyze and has little effect on the resulting 
luminosity function estimates. 

We now investigate explicitly the degree to which the 
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Figure 8. The points show the distribution of estimated cen- 
tral surface brightness, /io, and absolute magnitude, Mj, for our 
J< 14.45 (Kron) sample. The contours show visibility theory es- 
timates of Vmax as a function of /io and Mj. The contours are 
labelled by their Vmax values in units of (Mpc/h) 3 . 

completeness of the 2MASS catalogue depends on surface 
brightness by estimating Vmax as a function of both abso- 
lute magnitude and surface brightness. This is an impor- 
tant issue: if the catalogue is missing low-surface brightness 
galaxies our estimates of the luminosity function will be bi- 
ased. The approach we have taken follows that developed in 
Cross et al. (2001) for the 2dFGRS. We estimate an effec- 
tive central surface brightness, /iq, for each observed galaxy 
assuming an exponential light distribution, that the Kron 
magnitudes are total and that the Kron radii are exactly 
five exponential scale- lengths. This is then corrected to red- 
shift z = using 

/io =Ho - 101og(l + «)-A;(z)-e(z) (2.1) 

to account for redshift dimming and k+e corrections (c.f 
Section 3). Note that in the 2MASS catalogue, galaxies with 
estimated Kron radii less than 7 arcsec, have their Kron radii 
set to 7 arcsec. This will lead us to underestimate the central 
surface brightnesses of these galaxies, but this will only affect 
high surface brightness objects and will not affect whether a 
galaxy can or cannot be seen. The distribution in the Mj-/xo 
plane of our Kron J-band selected sample is shown by the 
points in Fig. 8. 

Cross et al. (2001) use two different methods to esti- 
mate the value of V ma x associated with each position in 
this plane. The first method uses the visibility theory of 
Phillipps, Davies & Disney (1990). We model the selection 
characteristics of the 2MASS extended source catalogue by 
a set of thresholds. The values appropriate in the J-band are 
a minimum isophotal diameter of 8.5 arcsec at an isophote 
of 20.5 mag arcsec -2 , and an isophotal magnitude limit of 
J<14.7 at an isophote of 21.0 mag arcsec -2 (Jarrett et al. 
2000). In addition, we impose the limits in the Kron magni- 
tude of 11<J<14.45 that define the sample we analyze. We 
then calculate for each point on the Mj-fio plane the red- 
shift at which a such a galaxy will drop below one or other 
of these selection thresholds and hence compute a value of 
Vmax- The results of this procedure are shown by the con- 
tours of constant V ma x plotted in Fig. 8. Note that these 
estimates of Vmax are only approximate since we have made 
the crude assumption that all the galaxies are circular expo- 




Mj-5 log h 

Figure 9. Three 1/Vmax estimates of the Kron J-band luminos- 
ity function. The data points with errorbars show the estimate 
based on assuming that Vmax depends only upon absolute mag- 
nitude and ignoring any possible surface brightness dependence. 
The dotted line and heavy solid line show the estimates in which 
the surface brightness dependence of Vmax is derived from visibil- 
ity theory and from the empirical method of Cross et al. (2000) 
respectively. 

nential disks. In addition, the diameter and isophotal limits 
are only approximate and vary with observing conditions. 

The second method developed by Cross et al. (2001) 
consists of making an empirical estimate of Vmax in bins in 
the Mj-fio plane. They look at the distribution of observed 
redshifts in a given bin and adopt the 90 th percentile of this 
distribution to define « max and hence V max . It is more robust 
to use the 90 th percentile rather than the 100 th percentile 
and the effect of this choice can easily be compensated for 
when estimating the luminosity function (Cross et al. 2001). 
Note that in our application to the 2MASS data we do not 
apply corrections for incompleteness or the effects of clus- 
tering. The result of this procedure is to confirm that for 
the populated bins, the V m ax values given by the visibility 
theory are a good description of the data. 

In Fig. 8 we see that the distribution of galaxies in the 
Mj-no plane is well separated from the low surface bright- 
ness limit of approximately 20.5 mag arcsec -2 where the 
Knax contours indicate that the survey has very little sen- 
sitivity. Thus, there is no evidence that low-surface bright- 
ness galaxies are missing from the 2MASS catalogue. Fur- 
thermore, in the region occupied by the observed data, the 
Knax contours are close to vertical indicating that there is 
little dependence of Kn ax on surface brightness. The way in 
which the V/ Vmax distribution is modified by including this 
estimate of the surface brightness dependence is shown by 
the dotted histogram in the top-left panel of Fig. 7. Its ef- 
fect is to increase the mean V/V ma x slightly, suggesting that 
this estimate perhaps overcorrects for the effect of surface 
brightness selection. Even so, the change in the estimated 
luminosity function is negligible as confirmed by the three 
estimates of the Kron J-band luminosity function shown in 
Fig. 9. These are all simple 1/V m ax estimates, but with V ma x 
computed either ignoring surface brightness effects or using 
one of the two methods described above. These luminosity 
functions differ negligibly, indicating that no bias is intro- 
duced by ignoring surface brightness selection effects. 
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Figure 10. The redshift distribution of the Kg<13.2 (Kron) sam- 
ple selected from the matched 2MASS-2dFGRS catalogue. The 
smooth curve is the model prediction based on the SWML es- 
timate of the Kg-band luminosity function (c.f. Section 4). The 
model prediction is very insensitive to the assumed k+e correction 
and cosmology. 



3 MODELLING THE STELLAR 
POPULATIONS 

The primary aim of this paper is to determine the present- 
day J and Ks-band luminosity functions and also the stellar 
mass function of galaxies. Since the 2MASS survey spans a 
range of redshift (see Fig. 10), we must correct for both the 
redshifting of the filter bandpass (k-correction) and for the 
effects of galaxy evolution (e-correction) . In practice, the k 
and e-corrections at these wavelengths are both small and 
uncertainties in them have little effect on the estimated lu- 
minosity functions. This is because these infrared bands are 
not dominated by young stars and also because the 2MASS 
survey does not probe a large range of redshift. We have cho- 
sen to derive individual k and e-corrections for each galaxy 
using the stellar population synthesis models of Bruzual & 
Chariot (1993; 2001 in preparation). We have taken this ap- 
proach not because such detailed modelling is necessary to 
derive robust luminosity functions, but because it enables 
us to explore the uncertainties in the derived galaxy stellar 
mass functions, which are, in fact, dominated by uncertain- 
ties in the properties of the stellar populations. 

The latest models of Bruzual & Chariot (2001 in prepa- 
ration) provide, for a variety of different stellar initial mass 
functions (IMFs), the spectral energy distribution (SED), 
lx(t,Z), of a single population of stars formed at the same 
time with a single metallicity, as a function of both age, t, 
and metallicity, Z. We convolve these with an assumed star 
formation history, ip(t'), to compute the time-evolving SED 
of the model galaxy, 



L\{t) = / lx(t-t',Z) dt'. 



(3-1) 



We take account of the effect of dust extinction on the SEDs 
using the Ferrara et al. (1999) extinction model normalized 
so that the V-band central face-on optical depth of the Milk- 
Way is 10. This value corresponds to the mean optical depth 
of Li, galaxies in the model of Cole et al. (2000) which em- 
ploys the same model of dust extinction. We assume a typi- 
cal inclination angle of 60 degrees which yields a net atten- 
uation factor of 0.53 in the V-band and 0.78 in the J-band. 
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Figure 11. The points in the upper two panels show the observed 
distributions of J-Kg and bj-Kgcolours as a function of redshift 
for our matched 2MASS-2dFGRS catalogue with z < 0.2 and 
J< 14.45 (Kron). Overlaid on these points are some examples of 
model tracks. The solid curves are for solar metallicity, Z = 0.02, 
and the dashed curves for Z = 0.004. Within each set, the tracks 
show different choices of the star formation time scale, r. The grid 
of values we use has r =1, 3, 5, 10 and 50 Gyr. Shorter values of r 
lead to older stellar populations and redder colours. The bottom 
panel shows the k+e corrections in the J-band for these same sets 
of tracks. 



By varying the assumed metallicity, Z, and star formation 
history, we build up a two-dimensional grid of models. Then, 
for each of these models, we extract tracks of bj— Ks and 
J— Ks colours and stellar mass-to-light ratio as a function of 
redshift. 

Our standard set of tracks assumes a cosmological 
model with fio = 0.3, Ao = 0.7, Hubble constant Ho = 70 
km s _1 Mpc -1 , and star formation histories with an expo- 
nential form, ip(t) oc exp(— [t(z) — t(zf)]/r). Here, t(z) is the 
age of the Universe at redshift z and the galaxy is assumed 
to start forming stars at Zf = 20 . For these tracks, we adopt 
the Kennicutt IMF (Kennicutt 1983) and include the dust 
extinction model. The individual tracks are labelled by a 
metallicity, Z, which varies from Z = 0.0001 to Z = 0.05 and 
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a star formation timescale, r, which varies from r = 1 Gyr to 
t — 50 Gyr. Examples of these tracks are shown in Fig. 11, 
along with the observed redshifts and colours of the 2MASS 
galaxies. We can see that the infrared J— Ks colour depends 
mainly on metallicity while the bj— Ks colour depends both 
on metallicity and star formation timescale. Thus, the use 
of both colours allows a unique track to be selected. Note 
from the bottom panel that, for all the tracks, the k+e cor- 
rection is always small for the range of redshift spanned by 
our data. 

We can gauge how robust our results are by varying the 
assumptions of our model. In particular, we vary the IMF, 
the dust extinction and cosmological models, and include or 
exclude the evolutionary contribution to the k+e correction. 
Also, we consider power-law star formation histories, ip(t) oc 
[t(z) /t(z{)]~"' , as an alternative to the exponential model. 
The results are discussed at beginning of Section 5. 

The procedure for computing the individual galaxy k+c 
corrections is straightforward. At the measured redshift of a 
galaxy, we find the model whose bj— Ks and J— Ks colours 
most closely match that of the observed galaxy. Having se- 
lected the model we then follow it to z = to predict the 
galaxy's present-day J and Ks-band luminosities and also 
its total stellar mass. We also use the model track to follow 
its k+e correction to higher redshift in order to compute 
Zmax, the maximum redshift at which this galaxy would have 
passed the selection criteria for inclusion into the analysis 
sample. 



4 LUMINOSITY FUNCTION ESTIMATION 

We use both the simple 1/V ma , x method and standard max- 
imum likelihood methods to estimate luminosity functions. 
We present Schechter function fits computed using the STY 
method (Sandage, Tammann & Yahil 1978) and also non- 
parametric estimates using the stepwise maximum like- 
lihood method (SWML) of Efstathiou, Ellis & Peterson 
(1988). Our implementation of each of these methods is de- 
scribed and tested in Norberg et al. (2001). The advantage 
of the maximum likelihood methods is that they are not 
affected by galaxy clustering (provided that the galaxy lu- 
minosity function is independent of galaxy density). By con- 
trast the 1/Vmax method, which makes no assumption about 
the dependence of the luminosity function with density, is 
subject to biases produced by density fluctuations. 

The two maximum likelihood methods determine the 
shape of the luminosity function, but not its overall normal- 
ization. We have chosen to normalize the luminosity func- 
tions by matching the galaxy number counts of Jarrett et 
al. (in preparation). These were obtained from a 184 deg 2 
area selected to have low stellar density and in which all 
the galaxy classifications have all been confirmed by eye. 
The counts are reproduced in Fig. 6. By using the same 
7 arcsec aperture magnitudes as Jarrett et al. (in prepara- 
tion) and scaling the galaxy counts in our redshift survey, 
we deduce that the effective area of our redshift catalogue is 
619 ± 25 deg 2 . Note that normalizing in this way by-passes 
the problem of whether or not some fraction of the missed 
2MASS objects are stars. Fig. 6 also shows the Kron and 
extrapolated magnitude J and Ks counts of the 2MASS- 
2dFGRS redshift survey. In the lower panel, these counts are 



seen to be in agreement with the published K-band counts 
of Gardner et al. (1996) and Glazebrook et al. (1994). 

We also checked the normalization using the following 
independent estimate of the effective solid angle of the red- 
shift survey. For galaxies in the 2dFGRS parent catalogue 
brighter than bj< Bn m i t , we computed the fraction that 
have both measured redshifts and match a 2MASS galaxy. 
For a faint Biimit this fraction is small as the 2dFGRS cat- 
alogue is much deeper than the 2MASS catalogue, but as 
-Biimit is made brighter, the fraction asymptotes to the frac- 
tion of the area of the 2dFGRS parent catalogue covered by 
the joint 2MASS-2dFGRS redshift survey. By this method 
we estimate that the effective area of our redshift catalogue 
is 642 ± 22 deg , which is in good agreement with the esti- 
mate from the counts of Jarrett et al. . 

It should be noted that for neither of these estimates 
of the effective survey area do the quoted uncertainties take 
account of variations in the number counts due to large scale 
structure. To estimate the expected variation in the galaxy 
number counts within the combined 2MASS-2dFGRS sur- 
vey due to large scale structure we constructed an ensemble 
of mock catalogues from the ACDM Hubble volume simula- 
tion of the VIRGO consortium (Evrard 1998; Evrard et al. 
in preparation; http://www.physics.lsa.umich.edu/hubble- 
volume). Mock 2dFGRS catalogues 

constructed from the VIRGO Hubble Volume simulations 
(Baugh et al. in preparation) can be found at http: //star- 
www. dur.ac.uk/ "cole/mocks/hubble. html . We simply took 
these catalogues and sampled them to the depth of 2MASS 
over a solid angle of 619 deg 2 . To this magnitude limit we 
found an rms variation in the number of galaxies of 15%. We 
took this to be a realistic estimate of the uncertainty in the 
2MASS number counts and propagated this error through 
when computing the error on the normalization of the lumi- 
nosity function. 



5 RESULTS 

5.1 Luminosity Functions 

Fig. 12 shows SWML estimates of the Kron J and Ks lumi- 
nosity functions. The points with errorbars show results for 
our default choice of k+e corrections, namely those obtained 
for an fi = 0.3, A = 0.7, H = 70 km s _1 Mpc -1 cosmol- 
ogy with a Kennicutt IMF and including dust extinction. 
The figure also illustrates that the luminosity functions are 
very robust to varying this set of assumptions. The various 
curves in each plot are estimates made neglecting dust ex- 
tinction and/or switching to a Salpeter IMF and/or chang- 
ing the Hubble constant to Ho = 50 km s _1 Mpc _1 and/or 
adopting power-law star formation histories and/or making 
a k-correction but no evolution correction. The systematic 
shifts caused by varying these assumptions are all compara- 
ble with or smaller than the statistical errors. The biggest 
shift results from applying or neglecting the evolutionary 
correction. In terms of the characteristic luminosity in the 
STY Schechter function fit, the estimates which include evo- 
lutionary corrections are 0.05 to 0.1 magnitudes fainter than 
those that only include k-corrections (see Table 1). 

In Fig. 13 we compare 1/V m ax and SWML Kron lumi- 
nosity function estimates (for our default choice of k+e cor- 
rections) with STY Schechter function estimates. In general, 
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Figure 12. SWML estimates of the Kron magnitude J (left) and Ks-band (right) luminosity functions (points with error bars). Our 
default model of k+e corrections (Kennicutt IMF and standard dust extinction) is adopted. The set of curves on each plot shows the 
effects of neglecting dust extinction and/or switching to a Salpcter IMF and/or changing the Hubble constant to Ho = 50 km s _1 Mpc _1 
and/or adopting power-law star formation histories and/or making a k-correction but no evolution correction. 



Table 1. The dependence of the J and Ks-band Schechter function parameters on cosmological parameters and evolutionary corrections. 
The parameters refer to STY estimates of the luminosity function for the 2MASS Kron magnitudes and are derived using k or k+c 
corrections based on model tracks that include dust extinction and assume the Kennicutt IMF. To convert to total magnitudes we 
estimate that the M* values should be brightened by between AKron + ^^j° n ~ A conv = 0.08 and 0.15 magnitudes. Note that the 
statistical errors we quote for the <I>* values include the significant contribution that we estimate is induced by large scale structure. 
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the luminosity functions are well fit by Schechter functions, 
but there is marginal evidence for an excess of very luminous 
galaxies over that expected from the fitted Schechter func- 
tions. We tabulate the SWML estimates in Table 2. Integrat- 
ing over the luminosity function gives luminosity densities 
in the J and K s -bands of p,/ = (2.75 + 0.41) x 10 8 /iL Q Mpc" 3 
and pK s = (5.74 + 0.86) x lO 8 /iL Mpc~ 3 respectively, where 
we have adopted Mf = 3.73 and M® = 3.39 (Allen 1973; 
Johnson 1966). In this analysis, we have not taken account 
of the systematic and random measurement errors in the 
galaxy magnitudes. In the case of the STY estimate, the 
random measurement errors can be accounted for by fitting 
a Schechter function which has been convolved with the dis- 
tribution of magnitude errors. However, for the Kron magni- 
tudes, the rms measurement error is only 0.1 magnitudes, as 
indicated by the comparison in the top right hand panel of 
Fig 5, and such a convolution has only a small effect on the 
resulting Schechter function parameters. We find that the 
only parameter that is affected is M* which becomes fainter 



by just A C onv = 0.02 magnitudes. The comparison to the 
Loveday (2000) data also indicates a systematic error in the 
2MASS Kron magnitudes of A^ r ° n = 0.061 + 0.031. Combin- 
ing these two systematic errors results in a net brightening of 
M* by AS on - A conv = 0.041 ± 0.031 magnitudes. As this 
net systematic error is both small and uncertain we have 
chosen not to apply a correction to our quoted Kron mag- 
nitude luminosity function parameters. We recall also that 
to convert from Kron to total magnitudes requires bright- 
ening M* by between AKron = 0.044 and 0.11 depending on 
whether the luminosity profile of a typical galaxy is fit well 
by an exponential or r 1/,4 -law. 

Fig. 14 shows the SWML and STY luminosity func- 
tion estimates for samples defined by the 2MASS extrapo- 
lated, rather than Kron, magnitudes. With this definition 
of magnitude, the luminosity functions differ significantly 
from those estimated using Kron magnitudes. In particular, 
the characteristic luminosities are 0.34 and 0.28 magnitudes 
brighter in J and Ks respectively. Most of this difference 
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Figure 13. SWML estimates of the Kron magnitude J (left) and Kg-band (right) luminosity functions (data points with errorbars) and 
STY Schechter function estimates (lines). The parameter values and error estimates of the Schechter functions are given in the legends. 
The errorbars without data points show 1/V max estimates of the luminosity functions. For clarity these have been displaced to the left 
by 0.1 magnitudes. 




Figure 14. SWML estimates of the extrapolated magnitude J (left) and Kg-band (right) luminosity functions (data points with errorbars) 
and STY Schechter function estimates (lines). The parameter values and error estimates of the Schechter functions are given in the legends. 
The errorbars without data points show 1/V ma x estimates of the luminosity functions. For clarity these have been displaced to the left 
by 0.1 magnitudes. 



is directly related to the systematic offset in the J-band 
Kron and extrapolated magnitudes, which can be seen in 
either the middle panel of Fig. 4 or the right hand panels 
of Fig. 5 to be approximately 0.23 magnitudes. Note that 
even in the Ks-band, it is this J-band offset that is relevant 
as the Ks-band magnitudes we use are derived from the J- 
band values using the measured aperture colours. We have 
argued in Section 2.3 that this offset is caused by the J-band 



Kron 2MASS magnitudes being fainter than true total mag- 
nitudes by between A^ r ° n + Apcron = 0.1 and 0.17 and the 
extrapolated magnitudes being systematically too bright by 
_ A oxtra P . _ A Kron = 0.05 to 0.11 magnitudes. Subtracting 
this 0.23 magnitude offset results in Kron and extrapolated 
luminosity functions that differ in M* by only 0.11 magni- 
tudes. In the Ks-band, the faint end slope of the best-fit 
Schechter function is significantly steeper in the extrapo- 
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Table 2. The SWML J and Kg-band luminosity functions for 
Kron magnitudes as plotted in Fig. 13. The units of both <j> and 
its uncertainty A<f> are number per h~ 3 Mpc 3 per magnitude. 
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lated magnitude case, but note that this function is not a 
good description of the faint end of the luminosity function 
since the SWML and l/Knax estimates lie systematically 
below it. The Schechter function fit is constrained mainly 
around M* and in this case the \ 2 value indicates it is not 
a good fit overall. 

The residual differences between the Kron and extrapo- 
lated magnitude luminosity functions arise from the scatter 
in the relation between extrapolated and Kron magnitudes. 
If this scatter is dominated by measurement error, then these 
differences represent small biases, which are largest for the 
less robust, extrapolated magnitudes. However, it is possi- 
ble that the scatter is due to genuine variations in galaxy 
morphology and light profiles. To assess which of these al- 
ternatives is correct requires independent deep photome- 
try of a sample of 2MASS galaxies to quantify the accu- 
racy of the extrapolated magnitudes. However, we note that 
the VyV m ax distributions for the extrapolated magnitudes 
shown in Fig. 7 have mean (V/Vmax) values significantly 
greater than 0.5, which is probably an indication that the 
extrapolated magnitudes are not robust. Thus, overall we 
favour adopting Kron magnitudes, noting the small offset of 
Af*° n + A Kr on - A C onv = 0.08 to 0.15 required to convert 



to total magnitudes and correct for the convolving effect of 
measurement errors. 

The parameters of the STY Schechter function fits 
shown in Fig. 13 are listed in the first row of Table 1. 
The subsequent rows illustrate how the best-fit parameters 
change when the cosmological model is varied and the evo- 
lutionary correction is included or excluded. The M* values 
are approximately 0.14 magnitudes fainter for the fio = 1 
case than for our standard flo = 0.3, Ao = 0.7 cosmology. 
This shift is largely due to the difference in distance mod- 
uli between the two cosmologies at the median redshift of 
the survey. This, and the difference in the volume-redshift 
relation, cause <j>* to change in order to preserve the same 
galaxy number counts. 

Fig. 15 compares our estimates of the Ks-band luminos- 
ity function for our standard fio =0.3, Ao =0.7 cosmology, 
with the estimates of Mobasher et al. (1993), Glazebrook et 
al. (1995), Gardner et al. (1997), Szokoly et al. (1998), Love- 
day (2000) and Kochanek et al. (2001). In general, these au- 
thors assumed different cosmological models when analysing 
their data. We have therefore modified the estimates from 
each survey. First, we apply a shift in magnitude reflecting 
the difference in distance moduli, at the median redshift, 
between the assumed cosmological model and our standard 
fio = 0.3, Ao = 0.7 model. We then apply a shift in number 
density so as to keep fixed the surface density of galaxies 
per square degree at the survey magnitude limit. In the case 
of the Kochanek et al. (2001) luminosity function we have 
shifted the data points brightwards by 0.05 magnitudes to 
account for the mean difference between the 2MASS isopho- 
tal magnitudes used by Kochanek et al. (2001) and the Kron 
magnitudes we have adopted. Schechter function parameters 
scaled and adjusted in this manner are given for each sur- 
vey in Table 3. Note that due to the correlations between 
the Schechter function parameters it is better to judge the 
agreement between the different estimates by reference to 
Fig. 15 rather than by the parameter values in Table 3. Our 
new estimate of the Ks-band luminosity function is in excel- 
lent agreement with the independent estimates and has the 
smallest statistical errors at all magnitudes brighter than 
Mk s — 5 log h = —22. For very faint magnitudes, from —20 
to —16, the sparsely sampled survey of Loveday (2000) has 
smaller statistical errors. Note that many previous analyses 
of the K-band luminosity function ignored the contribution 
of large scale structure to the error in and so the errors 
in Table 3 are likely to be underestimated. 

Also shown on the lower panel of Fig. 15 is an esti- 
mate of the K-band luminosity function inferred from the 
Sloan Digital Sky Survey (SDSS) near infrared, z*-band lu- 
minosity function of Blanton et al. (2001). To convert from 
z* (AB system) to standard K we have simply subtracted 
2.12 magnitudes from the SDSS z* magnitudes. This off- 
set consists of a contribution of 0.51 magnitudes to convert 
from AB magnitudes to the standard Vega system, and a 
mean z*-Ks colour of 1.61, which we find is typical of the 
model spectra discussed in Section 3 that match our ob- 
served b.j-Ks colours. As has been noted by Wright (2001) 
the luminosity function inferred from the SDSS data is off- 
set compared to our estimate. One suggestion put forward 
by Wright (2001) is that the 2MASS magnitudes could be 
systematically too faint. The systematic error would have 
to amount to 0.5 magnitudes to reconcile the luminosity 
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Figure 15. Comparison of various estimates of the K-band luminosity function. In both panels the solid line shows our SWML estimate 
of the Kg-band luminosity function for Kron magnitudes. The symbols and errorbars in the top panel show the estimates of Mobasher 
et al. (1993), Glazebrook et al. (1995), Gardner et al. (1997) and Szokoly et al. (1998) as indicated in the legend. We have shifted the 
data of Glazebrook et al. (1995) brightward by 0.3 magnitudes and the data of Mobasher et al. (1993) faintward by 0.22 magnitudes 
as advocated by Glazebrook et al. (1995) to make aperture corrections and consistent k-corrections respectively. In the lower panel the 
symbols and errorbars show the recent estimates of Loveday (2000) and Kochanek et al. (2000). The estimate of Kochanek et al. has 
been shifted brightward by 0.05 magnitudes to account for the difference between isophotal and Kron magnitudes. The dotted line shows 
a Schechter function estimate of the K-band luminosity function inferred from SDSS z*-band luminosity function of Blanton et al. (2000) 
(see text for details). 
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Table 3. Schechtcr function fits to K-band luminosity functions. Where necessary, the values quoted have been converted from the 
cosmological model assumed in the original work to the do = 0.3, Ao = 0.7 assumed here. In addition, we have shifted the M£ 
of Kochanek ct al. (2000) brightward by 0.05 magnitudes corresponding to the mean difference between 2MASS Kron and isophotal 
magnitudes in the Kochanek ct al. sample. We have also shifted M£ of Glazebrook et al. (1995) brightward by 0.3 magnitudes and that 
of Mobasher et al. (1993) faintward by 0.22 magnitudes as advocated by Glazebrook et al. (1995) to make aperture corrections and 
consistent k-corrcctions respectively. 
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Figure 16. The distribution of rest frame bj— Ks (left) and J— Kg (right) colours in three bins of Kg absolute magnitude, computed 
using our default set of k±e corrections. 



density inferred from the SDSS data with that which we 
infer from the 2MASS-2dFGRS catalogue. Such an error is 
comprehensively excluded by the very small offset that was 
found in Section 2.3 between the 2MASS Kron magnitudes 
and the data of Loveday (2000). Also, a direct galaxy-by- 
galaxy comparison of the z*-J and z*-Ks colours computed 
using the SDSS Petrosian and 2MASS Kron magnitudes pro- 
duced galaxy colours in good accord with expectations based 
on model spectra (Ivezic, Blanton and Loveday private com- 



munication). Finally, we note that a good match to our esti- 
mate of the Ks-band luminosity function cannot be achieved 
by simply moving the SDDS curve in Fig. 15 horizontally. 
If slid by 0.5 magnitudes to match the luminosity density 
then it falls well below our estimate at bright magnitudes. 
However if the SDSS curve is moved vertically, by a factor of 
1.6, then the two estimates come into reasonable agreement 
at all magnitudes. Thus, the most likely explanation of the 
difference between the SDSS and 2MASS-2dFGRS luminos- 
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ity functions is the uncertainty in the overall normalization 
induced by large-scale density fluctuations. It is to be hoped 
that as the sky coverage of the SDSS and 2MASS-2dFGRS 
surveys increases this discrepancy will be reduced. 



5.2 Colour Distributions 

Since our combined 2MASS-2dFGRS catalogue includes bj- 
band and infrared magnitudes, it is also possible to es- 
timate the bj-band optical luminosity function and the 
optical/infrared bivariate luminosity function. We do not 
present the bj-band optical luminosity function here as esti- 
mates from the 2dFGRS are discussed in detail in Norberg 
et al. (2001) and decomposed into luminosity functions of 
different spectral types in Folkes et al. (1999) and Madg- 
wick et al. (2001). Instead, we present the bivariate bj/Ks 
and J/Ks luminosity functions in Fig. 16, in the form of 
the rest-frame bj— Ks and J— Ks colour distributions, split 
by Ks-band absolute magnitude. Results are shown for just 
our default set of k+e corrections, but the colour distribu- 
tions are extremely insensitive to this choice and to whether 
evolutionary corrections are ignored or included. The shape 
of bj— Ks colour distribution varies systematically with Ks- 
band luminosity. At fainter magnitudes there is an increas- 
ingly large population of bluer, star-forming galaxies. The 
star formation rate has less effect on the infrared J— Ks 
colours. Here, the shape of the J— Ks colour distribution 
varies little with luminosity, but the position of the peak 
moves gradually redder with increasing luminosity. Colour 
distributions such as these are sensitive to both the distribu- 
tion of stellar age and the metallicity, and therefore provide 
important constraints on models of galaxy formation (for 
example, see Cole et al. 2000). 

5.3 Spectral Type Distribution 

Another interesting issue that we can address with our data 
is the distribution of spectral types in the 2MASS catalogue. 
For this, we make use of the spectral information in the 2dF 
galaxies extracted by a principal component analysis (Folkes 
et al. 1999). Specifically, we use the new continuous variable 
introduced in Madgwick, Lahav & Taylor (2001) which is 
defined by a linear combination of the first 2 principal com- 
ponent projections, r\ = 0.44pci — pc2- This variable was 
chosen to be robust to instrumental uncertainties whilst, at 
the same time, preserving physical information about the 
galaxy. The dominant influence on the r\ parameter is the 
relative strength of absorption and emission lines (?7<0 im- 
plies less than average emission-line strength while T7>0 im- 
plies stronger than average emission-line strength). A more 
detailed description will be presented in Madgwick et al. (in 
preparation). 

We can now gain insight into the population mix of our 
2MASS sample by simply creating a histogram of the r\ val- 
ues for the corresponding 2MASS-2dFGRS matched galax- 
ies with J< 14.45 (Kron). We plot this in Fig. 17 where we 
also show data for the entire 2dFGRS sample as comparison. 
Also shown in Fig. 17 (bottom panel) is the morphology-?7 
relation derived from a sample of galaxy spectra from the 
Kennicutt Atlas (Kennicutt 1992). 

It can be clearly seen from Fig. 17 that the predominant 
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Figure 17. The distribution of the spectral type parameter, 77, 
in the full 2dFGRS and our matched 2MASS catalogue (upper 
panel). The lower panel uses a sample of galaxies from Kennicutt 
(1992) to show how r\ is correlated with morphological type. 

population in the 2MASS sample is has r\<—2. By contrast- 
ing this with values of r\ obtained from the spectra of galax- 
ies with known morphological type (Kennicutt 1992) , we can 
see that this corresponds to galaxies of E/S0 morphologies. 
More precisely, the fraction of galaxies in our matched sam- 
ple with spectral types corresponding to E/S0 morphologies 
is 62% (compared with ~ 35% in the full 2dFGRS). Sa-Sb 
galaxies make up a further 22% and the remaining 16% are 
galaxies of later morphological types. 

5.4 The Galaxy Stellar Mass Function 

In contrast to optical light, near- infrared luminosities are 
relatively insensitive to the presence of young stars and can 
be more accurately related to the underlying stellar mass. 
Thus, with relatively few model assumptions, we can de- 
rive the distribution of galaxy stellar masses. The integral 
of this distribution is the total mass density in stars, which 
can be expressed in units of the critical density as O st ars ■ At- 
tempts to estimate this quantity date back many decades, 
but even recent estimates such as those by Persic & Salucci 
(1992), Gnedin & Ostriker (1992), Fukugita, Hogan & Pee- 
bles (1996) and Salucci & Persic (1999) have very large un- 
certainties because they are based on B-band light and re- 
quire uncertain B-band mass-to-light ratios. The much more 
accurate estimate that we provide here should prove very 
useful for a variety of purposes. 

To estimate the galaxy stellar mass function, we use the 
modelling of the stellar populations described in Section 3 to 
obtain estimates of the present luminosity and stellar mass- 
to- light ratio for each galaxy in the survey. This is done on a 
galaxy- by- galaxy basis as described in Section 3. The sample 
we analyze is defined by the 11<J< 14.45 (Kron) apparent 
magnitude limits. The stellar mass that we estimate for each 



© 0000 RAS, MNRAS 000, 000-000 



18 S. Cole et al. 




9 10 11 12 9 10 11 12 



log (M stars /h-2M ) log (M slars /h-2M ) 

Figure 18. SWML estimates of the stellar mass function (open symbols with error bars) and STY Schechter function estimates (lines). 
The parameter and error estimates of the Schechter function fits are given in the legends. The left-hand panel is for a Kennicutt IMF 
with recycled fraction R = 0.42 and the right-hand panel for a Salpeter IMF with R = 0.28. 



galaxy is the mass locked up in stars and stellar remnants. 
This differs from the time integral of the star formation rate 
because some of the material that goes into forming massive 
stars is returned to the interstellar medium via winds and 
supernovae. For a given IMF, this recycled fraction, R, can 
be estimated reasonably accurately from stellar evolution 
theory. Here, we adopt the values R = 0.42 and 0.28 for the 
Kennicutt (1983) and Salpeter (1955) IMFs respectively, as 
described in Section 5.2 of Cole et al. (2000) who made use of 
the models of Renzini & Voli (1981) and Woosley & Weaver 
(1995). Hence, the stellar masses we choose to estimate are 
(1—R) times the time integral of the star formation rate to 
the present day. Note that the IMFs we consider assume that 
only stars with mass greater than 0.1M© ever form and so 
we are not accounting for any mass that may be locked up 
in the form of brown dwarfs. 

Our results are presented in Fig. 18 which shows both 
SWML and Schechter function estimates of the present- 
day galaxy stellar mass function for two choices of IMF. 
The SWML estimates are tabulated in Table 4. Just as for 
the luminosity functions, the stellar mass function is quite 
well described by the Schechter functional form. Integrat- 
ing over these Schechter functions to determine the total 
stellar mass gives r2 stars /i = (1.4 ± 0.21) x 10" 3 for the 
Kennicutt IMF and fi 8 t arB /i = (2.6 ± 0.39) x 10~ 3 for the 
Salpeter IMF. Note that the integral converges rapidly at 
both limits and, in particular, the contribution to f2 s tars 
from objects with M st ars < 10 9 h -2 Mq is negligible. We 
find that these values vary by less than the quoted errors 
when we alter the assumed (k+e)-corrections by either ig- 
noring evolution, ignoring dust or changing Qo- Taken to- 
gether with our estimates of the Ks-band luminosity den- 
sity these, estimates imply mean stellar mass-to-light ra- 
tios of 0.73 Mq /Lq in the case of the Kennicutt IMF and 
I.32M0/L0 for the Salpeter IMF. If we apply the correc- 



Table 4. The SWML stellar mass functions as plotted in Fig. 18. 
The units of both <f> and its uncertainty A0 are number per h~ 3 
Mpc 3 per decade of mass. 
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Figure 19. Observational estimates of the star formation history 
of the universe. The points with error bars show estimates of the 
mean star formation rate per unit volume at various rcdshifts 
(see Steidel et al. 1999 and references therein). The solid sym- 
bols are the star formation rates implied if there is no absorption 
by dust. The open symbols show estimates corrected for dust 
absorption using a Calzctti (1999) extinction law with a mean 
E(B - V) = 0.15 (Steidel ct al. 1999). In both cases an Q = 0.3, 
Ao = 0.7 cosmology has been used to calculate the volume as a 
function of rcdshift and a Salpcter IMF to convert luminosity to 
star formation rate. The smooth curves are the fits we use when 
integrating over time to estimate the total mass density of stars 
formed by the present. 

tion we estimated in Section 2.3 to transform 2MASS Kron 
into total magnitudes, then these estimates and their un- 
certainties increase to fLtarsft = (1.6 ± 0.24) x 10~ 3 for the 
Kennicutt IMF and fi sta rs/i = (2.9 ± 0.43) x 10~ 3 for the 
Salpeter IMF. Both of these estimates are consistent with 
the value, s2 sta rs = (3.0 ± 1.0) x 10~ 3 , derived by Salucci & 
Persic (1999) but have fractional statistical errors which are 
several times smaller. With our method, the uncertainty in 
fistars is clearly dominated by the uncertainty in the IMF. 
For some purposes, it is not possible to improve upon this 
without a more precise knowledge of the true IMF - assum- 
ing there is a universal IMF. However, for other applications, 
such as modelling the star formation history of the universe, 
it is necessary to assume a specific IMF to convert the ob- 
servational tracers of star formation to star formation rates. 
Hence, in this case, it is the much smaller statistical errors 
that are relevant. 

It is interesting to compare our values with what 
is inferred by integrating the observational estimates of 
the mean star formation history of the universe. Fig. 19 
shows observational estimates for one particular choice 
of cosmology and IMF and illustrates how the rates are 
sensitive to the assumed dust extinction. By fitting a 
smooth curve through these estimates, we can calculate 
the mass of stars formed by the present day and how 
this depends on the IMF and assumed dust extinction. 
The upper smooth curve shown in Fig. 19 is of the 
form = (a + bz) / (I + (z/c) d )hM Q yT~ 1 Mpc~ 3 , where 
(a, b, c, d) = (0.0166, 0.1848, 1.9474, 2.6316). The data points 



Table 5. Estimates of the present-day mass in stars and stellar 
remnants obtained by integrating over observational estimates of 
the star formation history of the universe. We express this stellar 
mass density in terms of the critical density and give values of 
fistars^ 2 estimated for different assumed IMF's and dust correc- 
tions. All values are for an f2rj = 0.3, Ao = 0.7 cosmology and 
assume stellar populations of solar mctallicity. 

Dust Extinction Kennicutt IMF Salpeter IMF 

E(B-V)=0.05 0.80 x 10" 3 1.30 x 10" 3 

E(B-V)=0.10 1.17 x 10~ 3 1.86 x 10~ 3 

E(B-V)=0.15 1.63 x 10" 3 2.66 x 10" 3 



uncorrected for dust extinction are fit with (a, b, c, d) = 
(0.0,0.0798,1.658,3.105). As for our estimates above, we 
assume that no mass goes into forming brown dwarfs and 
multiply the star formation rate by 1— R, where R is the 
recycled fraction, so as to form an estimate of the mass 
locked up in stars. Values of f2 s tarsft estimated in this way 
are listed in Table 5. The values in this table are for an 
Slo = 0.3, Ao = 0.7 cosmology, but they are insensitive to 
this choice. They depend slightly on the assumed metallic- 
ity of the stellar population and would be 10% lower if half 
solar, rather than solar metallicity were assumed. Note that 
the Tatars values inferred from the star formation history 
of the universe scale differently with the assumed Hubble 
constant than those inferred above from the IR luminosity 
functions. For h = 0.7 our estimates from 2MASS become 
n stals h 2 = (1.12 ± 0.16) x 10~ 3 for the Kennicutt IMF and 
n st arsk 2 = (2.03 ± 0.30) x 10~ 3 for the Salpeter IMF. Com- 
parison with Table 5 shows that these values are consistent 
with those inferred from the cosmic star formation history 
only if the dust correction assumed in the latter is mod- 
est, E(B-V)« 0.1. This value is 50% smaller than the value 
preferred by Steidel et al. (1999). 



6 CONCLUSIONS 

The new generation of very large surveys currently underway 
make it possible to characterize the galaxy population with 
unprecedented accuracy. In this paper, we have combined 
two such large surveys, the infrared imaging 2MASS and the 
2dF Galaxy Redshift Survey^ to obtain a complete dataset 
which is more than an order of magnitude larger than previ- 
ous datasets used for statistical studies of the near-infrared 
properties of the local galaxy population. We have used this 
combined catalogue to derive the most precise estimates to 
date of the galaxy J and Ks-band luminosity functions and 
of the galaxy stellar mass function. 

Characterizing the near-infrared properties of galax- 
ies offers several advantages. Firstly, the near-infrared light 
is dominated by established, old stellar populations rather 
than by the recent star formation activity that dominates 
the blue light. Thus, the J and K-band luminosity functions 
reflect the integrated star formation history of a galaxy and, 

T A table containing the positions, 2MASS infra-red magni- 
tudes and 2dFGRS redshifts used in this paper and elec- 
tronic versions of Tables 2 and 4 are available at http://star- 
www.dur.ac.uk/-cole/2dFGRS-2MASS . 
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as a result, provide particularly important diagnostics of the 
processes of galaxy formation. For the same reason, the dis- 
tribution of stellar mass in galaxies -the galaxy stellar mass 
function- can be derived from the near-infrared luminosities 
in a relatively straightforward way, with only a weak model 
dependence. Finally, corrections for dust extinction as well 
as k-corrections are much smaller in the near-infrared than 
in the optical. 

Due to the size of our sample, our determination of 
the J- and Ks-band galaxy luminosity functions have, for 
the most part, smaller statistical errors than previous esti- 
mates. Furthermore, since our sample is infrared-selected, 
our estimates are free from any potential biases that might 
affect infrared luminosity functions derived from optically- 
selected samples. We find that the J- and Ks-band galaxy 
luminosity functions are fairly well described by Schechter 
functions, although there is some evidence for an excess of 
bright galaxies relative to the best-fit Schechter functional 
form. In general, the SWML estimates are a truer represen- 
tation of the luminosity functions. Our K-band estimates are 
in overall agreement with most previous determinations, but 
have smaller statistical errors. 

The exception is the K-band luminosity function in- 
ferred from the near infrared SDDS photometry (Blanton 
et al. 2001). The difference between the K-band luminosity 
function we infer from their data and our own estimate is 
too large to be explained by photometric differences. The 
difference between the two estimates is better described by 
a difference in overall number density of a factor of 1.6. 
A similar discrepancy is seen in the bj-band between the 
SDSS and 2dFGRS luminosity function estimates (see Nor- 
berg et al. 2001 and Blanton et al. 2001). The suspicion is 
that the uncertainty in the overall normalization of the lumi- 
nosity functions induced by large-scale structure within the 
large, but finite, survey volumes could be to blame. However, 
the errors that we quote for the 2dFGRS-2MASS luminosity 
functions already include an estimate of this sampling un- 
certainty as derived from realistic mock catalogues. A sim- 
ilar exercise for a catalogue with the same area and depth 
as that of Blanton et al. (2001) indicates that the required 
overdensity of a factor of 1.6 is unlikely. So probably there 
is more than one contributory factor at work and the hope 
is that these will be identified as the surveys progress. 

Using our J-band luminosity function, bj— Ks and 
J— Ks colours and simple galaxy evolutionary tracks, we 
have obtained the first estimate of the galactic stellar mass 
function derived directly from near-infrared data. We find 
that this mass function is also fairly well described by a 
Schechter form. An integral over the stellar mass function 
gives Tatars, the universal mass density locked up in lumi- 
nous stars and stellar remnants, expressed in terms of the 
critical density. n sta rs is a key component of the overall in- 
ventory of baryons in the universe. An accurate determi- 
nation of this quantity is essential for detailed comparisons 
with other quantities of cosmological interest such as the 
total baryonic mass, fibaryon, inferred from Big Bang nucle- 
osynthesis considerations (e.g. Buries & Tytler 1998), the 
cosmic star formation rate (e.g. Steidel et al. 1999), and the 
cosmic evolution of the gas content of the universe (Storrie- 
Lombardi & Wolfe 2000). The statistical uncertainty in 
our estimate of f2 s tars is about 15%, several times smaller 
than the best previous determination by Salucci & Persic 



(1999). In fact, the errors in f2 s tars are dominated by sys- 
tematic uncertainties associated with the choice of stellar 
IMF in the galaxy evolution model. A Kennicutt IMF gives 
n stals h = (1.6 ± 0.24) x 10~ 3 while a Salpeter IMF gives 
fWs/i = (2.9 ± 0.43) x 10~ 3 . For h = 0.7, these values cor- 
respond to less than 11% of the baryonic mass inferred from 
Big Bang Nucleosynthesis (Buries et al. 1999). Our values of 
fistars today are only consistent with recent determinations 
of the integrated cosmic star formation if the correction for 
dust extinction is modest. 
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